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Abstract —This paper presents a generalized energy storage 
system model for voltage and angle stability analysis. The 
proposed solution allows modeling most common energy storage 
technologies through a given set of linear differential alge¬ 
braic equations (DAEs). In particular, the paper considers, but 
is not limited to, compressed air, superconducting magnetic, 
electrochemical capacitor and battery energy storage devices. 
While able to cope with a variety of different technologies, the 
proposed generalized model proves to be accurate for angle and 
voltage stability analysis, as it includes a balanced, fundamental- 
frequency model of the voltage source converter (VSC) and 
the dynamics of the dc link. Regulators with inclusion of hard 
limits are also taken into account. The transient behavior of 
the generalized model is compared with detailed fundamental- 
frequency balanced models as well as commonly-used simplified 
models of energy storage devices. A comprehensive case study 
based on the WSCC 9-bus test system is presented and discussed. 

Index Terms —Energy storage system (ESS), transient stabil¬ 
ity, power system dynamic modeling, electrochemical capacitor 
energy storage (ECES), superconducting magnetic energy stor¬ 
age (SMES), compressed air energy storage (CAES), battery 
energy storage (BES). 

I. Introduction 

The capability of Energy Storage Systems (ESSs) to im¬ 
prove the transient behavior of power systems and to increase 
the competitiveness of non-dispatchable power production 
technologies, e.g., renewable energies, has led, in recent years, 
to a huge investment in research, prototyping and installation 
of storage devices as well as the development of a huge 
variety of energy storage technologies [l]-[3]. Among most 
relevant technologies that currently appear most promising 
there are Battery Energy Storage (BES), Compressed Air 
Energy Storage (CAES), Superconducting Magnetic Energy 
storage (SMES), Electrochemical Capacitor Energy Storage 
(ECES), and Elywheel Energy Storage (EES). 

Modeling ESSs is a complex and time-consuming task due 
to the number of different technologies that are currently 
available and that are expected to be developed in the future. 
While there are several studies aimed to define the economic 
viability and the effect on electricity markets of ESSs (see 
for example [2] and [4]), there is still no commonly-accepted 
simple yet accurate general model of ESSs for voltage and 
angle stability studies. This paper presents a generalized 
model of ESSs to simplify, without giving up accuracy, the 
simulation of different storage technologies. The objective is 
to provide a balanced, fundamental frequency model that can 
be defined through a reduced and fixed set of parameters and 
that can be readily implemented in power system simulators 
for voltage and angle stability analysis. 
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TABLE I: Examples of energy storage technologies. 


Types of 
Storable Energy 

Potential Var. 

Flow Var. 

Device 

Magnetic 

Magneto Motive 
Force 

Flux 

SMES 

Fluid 

Pressure 

Mass Flow 

CAES 

Electrostatic 

Electric Potential 

Electric Current 

ECES 

Electrochemical 

Electrochemical 

Potential 

Molar Flow 

Rate 

BES 

Rotational 

Angular Velocity 

Torque 

FES 


Studies that focus on ESS dynamic models aimed to tran¬ 
sient stability analysis abound in the literature. The following 
are relevant works. A general taxonomy for ESS dynamic 
analysis is presented in [5]. A detailed transient stability 
model of ECES can be found in [6], while the model of SMES 
is given in [7]. Studies of the dynamic performance of a SMES 
system coupled to a wind power plant throughout a voltage 
source converter (VSC) are presented in [8] and [9], whereas 
information regarding the installation of a large-scale SMES 
at the Center of Advanced Power Systems at Elorida State 
University can be found in [10]. In [11], the authors propose 
the model of a small-size CAES based on a polytropic process 
of air compression and expansion and in [12], [13] BES 
dynamic models are proposed and tested. Einally, transient 
stability models of EES can be found in [14] and [15]. 

Several simplified ESS models have been also proposed. 
Among these, we cite [16]-[20]. The main feature that the 
models proposed in the references above have in common is 
that the ESS dynamics are represented considering only active 
and reactive power controllers. However, the dynamics of the 
ESS itself are not preserved. These models are also generally 
loss-less as dc and VSC circuitry is neglected. 

In this paper, we use the balanced, fundamental frequency 
model of the VSC that is proposed in [21]-[25], which 
includes dc circuit and phase-locked loop (PEL) dynamics 
as well as an average dg-axis model of the converter and an 
equivalent model for switching losses. 

The model proposed is based on the observation that most 
ESSs connected to transmission and distribution grids share a 
common structure, i.e., are coupled to the ac network through 
a VSC device, present a dc-link and then include another 
converter (either a dc/dc or a ac/dc device) to connect the main 
energy storage device to the dc link. Moreover, all storage 
systems necessarily imply potential and flow quantities (see 
Table IT]), whose dynamics characterize the transient response 
of the ESS. In this paper, we show that, to properly capture 
the dynamic response of the ESS, it is important to preserve 
such dynamics along with those of the VSC converter and its 
controllers. Controller hard limits, whose relevance has been 
discussed in [20] and [26], are also considered. 

In summary, the paper provides the following contributions: 

1) A simple yet accurate generalized storage device model 
that consists of a given set of parameters and linear 
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differential algebraic equations (DAEs). This model is 
characterized by a given set of DAEs whose structure 
and dynamic order is independent from the technology 
used for the ESS. 

2) A balanced, fundamental frequency model of the typical 
devices that compose an ESS, namely, the VSC, the 
dc link and their regulators. This model is obtained 
based on a careful selection of well-assessed models 
for transient stability analysis. 

3) A comprehensive validation of the dynamic response 
of the proposed ESS model versus detailed transient 
stability models representing specific technologies as 
well as simplified models that include only active and 
reactive power controllers. 

The paper is organized as follows. Section |II] outlines 
the structure of a generic VSC-based ESS. The schemes of 
the regulators that are common to all ESSs, as well as the 
active and reactive controllers used for simplified ESS models 
are also presented in this section. Section [111] describes the 
hypotheses, the structure and the formulation of the proposed 
generalized model of ESS, as well as how this model is able to 
describe the behavior of a variety of ESS technologies. Section 
HVl compares the dynamic response of the proposed model 
with detailed and simplified models of ESSs for transient 
stability analysis. All simulations are based on the WSCC 
9-bus test system. Einally, section E draws conclusions and 
outlines future work directions. 
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Fig. 1: Scheme of the ESS connected to a grid, 
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Fig. 2: VSC scheme. 


II. ESS Controls and Simplified ESS Model 

This section presents the elements and controllers that 
regulate the dynamic response of the ESS. In this section, 
we consider exclusively the elements of the ESS model that 
are independent from the storage technology. We also present 
the simplified ESS model that is used for comparison in the 
case study presented in Section HVl 

A. Detailed ESS scheme for Transient Stability Analysis 

Eigure[T] shows the overall structure of an ESS connected to 
a grid through a VSC. The objective of the ESS is to control 
a measured quantity w, e.g., the frequency of the Center of 
Inertia (COI) or the power flowing through a transmission line. 
The elements shown in Eig. [T]are common to all VSC-based 
ESSs. The models of these elements, except for the storage 
device, are presented in the remainder of this section. 

Eigure |2] illustrates the VSC scheme [21]-[23]. The usual 
configuration includes a transformer, a bi-directional converter 
and a condenser. The transformer provides galvanic insulation, 
whereas the condenser maintains the voltage level at the dc 
side of the converter. The dc voltage is converted to ac voltage 
waveform via the use of power electronic converters that are 
controlled using appropriate control logic. 

Eigure 0 shows the block diagram representation of both 
the converter and the inner current control loop of the VSC. 
After transforming the three-phase equations to a rotating dq- 
frame, the dynamics of the ac side of the VSC depicted in 
Eig. |2] become; 


Inner Cunent Control Loop Converter 



Fig. 3: Block diagram of the inner current control and the converter 
in the dg-frame. 


where Rac+jLac is the aggregated impedance of the converter 
reactance and transformer; wt,d = m-dUdc/2; Ut,q = mqUdc/2; 
and Ldac is the frequency of the grid voltage, Fac. 

It is possible to express dc signals in a dg-reference frame, 
with the aim of designing suitable controllers, provided that 
the dg-reference frame rotates with the same speed as the grid 
voltage phasor, which may be achieved, for example, through 
the use of a phase-locked loop (PEL) [24]. The PEL is a 
control system that forces the angle of the rotating dg-frame, 
6dq, to track the angle dac- The PEL is typically composed 
of a phase detector (PD), a loop filter (EE), and a voltage 
controlled oscillator (VCO), and its scheme is depicted in 
Eig. a 
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Fig. 4: PLL scheme. 
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The grid voltage Va.c can, thus, be expressed in the dg-frame 
as follows; 

t^ac,d j'^ac,q — t^ac(cOs(dac ^dq) “t“ j sin(dac ^dq)^ (2) 

The controller iTi(s) of Fig. [3] is designed as a PI control, 
as follows [21]: 


PI control Storage Input Limiter 



l \ si/ac ,,, 

Ki{s) = - - - (3) 

sTi 

where Tj is the time-constant of the closed-loop step response. 

The power balance between the dc and the ac sides of the 
converter is imposed by: 


Pac “t” Ploss 2 


^(^dc) 

dt 


= 0 


(4) 


where pac = 


,d^a 


i); are the 

2 dt 


energy variations in the capacitor; and pioss = ^^adtc + 
Gsw(*dc)'fdc circuit and switching losses of the con¬ 
verter, respectively, with ^ -f and Gsw(*dc) 

is obtained from a given constant conductance, Gq, and the 
quadratic ratio of the actual current to the nominal one, as 
follows [25]: 


Gsw(*dc) — Go 


^dc 


Me 


(5) 


variable of the storage device, u, as follows: 


u = { 


E — E-n 


-Au + Mref 


rpthr _ TT . 

^min ^min 

if £^min <E < and Aw > 0 

E — E 

-^max A 

Am -f Mret 


(6) 


/? _ /^thr 

-^max -^max 


if <E< Fornax and Am < 0 


-thr 
ma: 

, M Otherwise 


where Mref is the value of u such that the storage device is 
disabled; Am = u — Mref; ^-min and E’max are the minimum 
and maximum storable energy in the ESS, respectively; and 
E^t^ and E^l^ are the given minimum and maximum energy 
thresholds, respectively. 

A dead-hand block is also included in Fig. |6l This block 
is designed to reduce the sensitivity of the Storage Control to 
small changes of the measured signal w with respect to the 
reference The purpose of the dead-hand is to reduce the 
number of charge/discharge operations of the ESS [9]. 


Finally, the outer controllers used to regulate the dc and ac 
voltages of the VSC are depicted in Fig. |5] 


PI control 



B. Simplified ESS scheme 

As discussed in SectionlH several simplified models of ESSs 
have been proposed [16]-[20]. For the sake of comparison, in 
this paper, we consider the control scheme proposed in [16] 
(see Fig.|2l). The input signal w is regulated through the active 
power while the voltage at the point of connection with the 
grid is regulated through the ESS reactive power. The physical 
behavior of the storage system is synthesized by the two lag 
blocks with time constants Tp^ESS and Tq ess- 



Fig. 5: Outer dc and ac voltage controllers. 


The charge/discharge process of the storage device is regu¬ 
lated by the storage control (see Fig. |6]l. The input signal of 
the control is the error between the actual value of a measured 
quantity of the system, say w, and a reference value (w’’®^). 
If w = w’'®^ the storage device is inactive and its energy is 
kept constant. For w ^ w’’®^, the storage device injects active 
power into the ac bus through the VSC (discharge process) or 
absorbs power from the ac bus (charge process). The scheme 
shown in Fig. |6] also includes a block referred to as Storage 
Input Limiter [26]. The purpose of this limiter is to smooth the 
transients that derive from energy saturations of the storage 
device. This block takes the actual value of the energy stored 
in the device, E, and defines the value of the controlled input 


PI control 



Fig. 7: Scheme of a simplified ESS model. 


III. Proposed Generalized Model oe Energy 
Storage Devices 

The generalized model proposed in this paper is based on 
the linearization and dynamic reduction of the set of equations 
that describes the storage device. Let assume that the dynamic 
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behavior of the ESS can be described by a set of nonlinear 
DAE as follows; 


Then, substituting (fTOl) into the first and third equations of 
(Ell, one has: 


ri(f) = 

= (7) 

where ^{t) is the state vector S R"); v{t) is the input 
vector {u G R™); ip{t) is the output vector (<p G R^); ip : 
R"+™ —>• R" are the differential equations; rj : R"+™ — 
RP are the output equations; and E is a diagonal matrix of 
dimensions n x n such that [27]: 

Tii = 1 if the i-th equation of ip is differential; 

Tii = 0 if the i-th equation of ip is algebraic. 

The first step towards the definition of the generalized 
model of the ESS is to linearize the system around the 
equilibrium point (^g, 93 q, zzo). Thus, the expression for a 
linear time-invariant dynamical system is obtained: 

ri(f) = A^{t) + Bzz(f) + Kj 
ipit) = C^{t) + -Du{t) + K^ (8) 

where A is the state matrix (dim [A] = n x n)-, B is 
the input matrix (dim[B] = n x m)', C is the output 
matrix (dim[C] = p x n); D is the feedthrough matrix 
(dim[D] = p X to); and Kj G R" and G R^ account for 
the values of the variables at the equilibrium point. Equation 
(O is written for u and <p, not the incremental values A^, 
Azz, and A<p. 

Then, the state vector ^ is split into two types of variables; 
the potential and flow variables related to the energy stored 
in the ESS, x (see Table HI; and all other state variables, z. 
Hence, ^ = \x'^, . Since all ESSs are coupled to the 

ac network through a VSC device, the dc voltage and current 
of the VSC, Udc and idc, are considered as an input and an 
output in dH, respectively. The input vector iz is composed 
of the output signals of the storage control, u, and the dc 
voltage of the VSC, Vdc- Hence, ix = [u, Udc]^. Einally, the 
output vector ip is the dc current of the VSC, *dc- Hence, 
V = [*dc] ■ Applying the notation above, equation ([^ is written 
as follows: 

= Axx^ “b B^^^W -f B^^^Udc “b 

TzZ = AzxX + AzzZ + B^uU + Bz„Udc + Kz (9) 
2dc — T CzZ -t- “b Da;Udc “b Ki 


where; 

A A 

■^xx -^xz 

A A 

■^zx -^zz 

[C, Cz]=C; [D„D„] = D; 

[r, FzY = r ; [K, Kzf = Kj ; K, = 

The dynamic order of (0 is reduced by assuming that 
the transient response of z are much faster than that of x 
or immaterial with respect to the overall dynamic behavior 
of the ESS. This assumption is based on the knowledge of 
detailed transient stability models and is duly verified through 
the simulations presented in Section IIVI By neglecting the 
dynamics of z (i.e., = 0), from the second equation of 

El, we obtain: 


Ba 

Bz 


= B 


BxX — {^Axx ^xz-^zz ^zx) ^ “b (Bxu ^xz-^zz Bzu) H 
“b (Ba;a, AxzA^z B^i;) Udc “b 
tdc = (Ca; — CzA^^^Aza;) X + (D„ — CzA^^^Bzu) U 

+ (D„ - CzA-iBz„) Udc + (Ka - CzA-iRz) (11) 

Rewriting in compact form the matrices in (fTTI . the pro¬ 
posed generalized model of energy storage devices can be 
written as follows: 

Ttc =Ax + B„u + B„Udc + Ka; 

tdc — Cx -f DatU “b E)a;Udc “b Ka (17) 

It is important to take into account the state of charge of the 
storage device (in particular, to impose energy limits to the 
controller shown in Eig. |6|. With this aim, the actual stored 
energy is given by; 

n 

7^ = EPi(^f‘-xf‘) (13) 

i=l 

where p;, Pi and Xi are the proportional coefficient, exponen¬ 
tial coefficient and reference potential value of each variable 
Xi, respectively. 

In summary, the proposed generalized ESS model is com¬ 
posed of the block diagram of the VSC represented in Eig. [2 
the controllers depicted in Eigs. |5]and|6j and (IT2l) . (fOl . 

Eor the sake of clarity and example, in the following 
subsections we deduce the generalized model for a variety 
of ESSs technologies. 

A. Electrochemical Capacitor Energy Storage Model 

Eigure [8] shows the scheme of an ECES connected to the 
VSC through a bidirectional dc/dc converter (buck-boost) [6]. 
In the buck operation mode, the energy is delivered from the 
storage device to the grid, while in the boost mode, the energy 
is stored in the ECES. 



The model that describes both the buck and the boost 
operation modes is as follows: 

ttsc— (^sc “b fzsct^sc) 

^SC 

^sc — "7 (t^sc tsc7(sc l^dcS) (14) 

-t-'sc 

Idc— Si^c 

where S is the logic state of the switches of the converter. 
Assuming average variable values, S is continuous and re¬ 
presents the duty cycle of the converter. 

The energy stored in the ECES is as follows: 

^ — 2^scl'sc + 


Z — ^zz {■^zx'^ ~b Bzull ~b Bzi;UUc ~b Kz) 


(10) 


(15) 
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In (fTSl) . the main term is that related to the capacitance Cgc, 
as expected. However, for completeness, we include also the 
energy stored in the inductance Lgc- 

Applying the notation of the proposed general model (fl^ 
m to (1141115b . ECES variables and parameters are; 


X = u 


SC ^scj ? 


: z = f 


= S: T = 


1 0 
0 1 


A = 

_ Gsc _ 1 

Csc C'sc 

; Bu = 

^ '^^dc,0 

T 

; Bu = 

[o ^1 


1 _ Rsc 

L/SC -i'SC , 


-^sc 




n T 


K, = 


{Gsc'^sc.O '^sc,o) 

Y, ( ^sc,0 “1“ ^sc^sc,0 O'S'o) 

L J-/SC ' ’ ’ ’ ' 

— [tsc,o]; — [0] 5 — [^dc,0 25^otsc,o]3 


; C = [0 5o]; 


P = 


-(7sc -As, 
2 2 


; /3 = [2 2]; X = [0 0] 


where Sq, Vsc,o, t^dc.o^ Ac,o and Zdc.o are the values of S, Vsc, 
Vdc, Ac and *dc at the equilibrium point, respectively. 


B. Superconducting Magnetic Energy Storage Model 

Figure 0 depicts the scheme of a SMES [7]. The SMES is 
connected in parallel to the VSC throughout a dc/dc converter 
(boost converter). The SMES stores magnetic energy and 
injects it into the network according to the duty cycle applied 
to the converter. 



Fig. 9: Scheme of the SMES and the boost converter. 



Atd.c md,t 

Fig. 10: Scheme of a Compressed Air Energy Storage. 


The air is modeled as an ideal gas. Therefore, the variation 
of the pressure inside the tank can be calculated as follows: 


tl2 


TTinVol^ 


(19) 


where 112 is the pressure inside the tank; p is the air density; 
R is the ideal gas constant; 02 is the temperature of the 
air inside the tank; is the molecular weight of air; Vol 
is the volume of the tank; and Q is the air flow through 
the compressor/turbine (all quantities are assumed in absolute 
value). 

Assuming that both compression and expansion of air are 
polytropic processes [11], the behavior of the compressor 
of the CAES is described by the following set of nonlinear 
equations (turbine equations are similar and are omitted): 


0 = 


7-1 


HiQ 


Ih 

Hi 


0 — Pei VmPm 
P 

0 — — — T 

0 = Tel — Pel 

„ Win — W2 

0 =-cr 

Win 

0 = Upfl — Ul2 


- 1 


-Pei 


( 20 ) 

( 21 ) 

( 22 ) 

(23) 

(24) 

(25) 

(26) 


The dynamics of the coil and the dc/dc converter are 
represented by the circuit equations: 


Ac 

0 = (1 - 2S)vdc - Vc 

idc= (1 - ‘2S)ic 


(16) 


where all voltages and currents are average values; Ac is the 
inductance of the superconducting coil; and S is the duty cycle 
of the dc/dc converter. 

The energy stored in the SMES is given by: 

E = (17) 

Therefore, using the notation of the proposed model: 


X = [A Vc] ; 2 = 0; u = S; T = 


1 0 
0 0 


p = 


Ilo 0 
2 


; /3 = [2 1]; x = [0 0] 


(18) 


C. Compressed Air Energy Storage Model 

The basic configuration of a CAES is represented in Fig. 
[TOl [11]. The air is injected into the tank by means of a 
compressor operated by an asynchronous motor, and extracted 
from it through a turbine driven by an asynchronous generator. 
Both the compressor and the turbine are connected in parallel 
to the dc link of the VSC through ac/dc converters. 


where 7 is the polytropic exponent of the air; Hi is the atmos¬ 
pheric pressure; Pef, Pm and Pei are the effective, mechanical 
and electrical powers of the compressor, respectively; Tm and 
Tel are the mechanical and electrical torques of the motor, 
respectively; Cl and u :2 are the mechanical and electrical rotor 
speed of the compressor, respectively; a is the motor slip; and 
77ni, rip, Win, H are motor parameters. 

For the sake of space, the equations of the electrical 
machines connected to the compressor and turbine are not 
shown in (i20l l- (i26l l. The model considered for the electric 
machines is a 5*^ order dg-model (see, for example. Chapter 
4 of [28]). Finally, the ac/dc converters are controlled in a 
similar way as the VSC described in Section HI! 

According to (i20l) and (i^Tb . the energy stored in the CAES 
can be expressed as follows: 


E = —^HiVol 
7- 1 



+ HCl^ 


(27) 


As in Subsection IIII-Al both terms related to the energy 
stored in the CAES have been considered, despite the fact that 
the term related to the mechanical rotor speed Cl is expected 
to be much smaller than the term related to the pressure 112. 

Taking into account that the compressor and the turbine 
have similar sizes, we assume that both compressor and 
turbine electrical machines have same parameters. Therefore, 
we model only one equivalent electrical machine capable of 
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working, depending on the sign of Q, as a compressor (Q > 0) 
or a turbine (Q < 0). 

Imposing the proposed generalized model to CAES equa¬ 
tions, we obtain: 


^ — [n2 11] , U — , tdc — tdc,c ^dc,t? 


2! = [Pef Pm Pel Pm Pel O’ 012 rhd TOq . . . (28) 


—Voi H 

7-1 


/3 


7-1 

7 



pi 0] 


where the dots in the vector 2 ; stand for electrical machine and 
VSC variables which are not explicitly defined in (I20l i- (l26l l. 
2 : contains 25 state and algebraic variables. 

Since we do not model the air valve, the output of the 
storage control, and thus the input to the generalized model 
u is the air flow, Q. By using the proposed model, the 
nonlinear set of DAEs which models the CAES ( (fT9l l- (l26l l. 
electrical machine equations of the compressor and turbine, 
and equations of the ac/dc converters) is reduced to a set of 
three linear DAEs (fTST i. 


D. Battery Energy Storage Model 

A commonly-used model to represent the dynamics of a 
rechargeable battery cell is the Shepherd model [13]: 


Qe = *b /3600 

*m — 

d m 

0 = Uoc t7p((^e5 Im) “f UeC Pi^b 


(29) 


Vb 


where Qe is the extracted capacity in Ah; is the battery 
current ib passed through a low-pass Alter with time constant 
Pm; Voc, Vp and Ve are the open-circuit, polarization and 
exponential voltages, respectively; /3e is the exponential zone 
time constant inverse; Pj is the internal battery resistance; and 
Ub is the battery voltage. 

The variation of the polarization voltage Vp with respect to 
im and Qe is given by: 


ttp {Qei ^m) — 


Pp^m ^pQe 

SOC 

Pp^m ^ ^pQe 


if 


ge + 0.1 SOC 


> 0 (discharge) 

if im <0 (charge) 
(30) 


where Pp and Kp are the polarization resistant and polariza¬ 
tion constant, respectively; and SOC is the state of charge of 
the battery which is defined as: 


SOC = 



= 1 - ge 


(31) 


The non-linearity of Vp implies that depending on the state 
of the battery (charge or discharge), two different sets of equa¬ 
tions can be obtained by applying the proposed generalized 
model. Therefore, the proposed model has to be able to switch 
from one set to another depending on the BES operation. 

The connection of the BES to the VSC is similar to the 
SMES (see Subsection IIII-Bb : 


0 = (1 - 2S)vdc - nsVh 

idc = - (1 - 2S)npii, (32) 

where S is the duty cycle of the converter; and rip and Ug 
are the number of parallel and series connected battery cells. 


respectively. The dc/dc converter connection used in this paper 
is based on [29]. Other configurations are also possible, but 
are not considered in this paper. 

Imposing the proposed generalized model to BES equa¬ 
tions, one has: 


— [Qe t^b] ; ^ — [^b Im t^p] 7 ^ — ^7 ^ — 


1 0 
0 0 


p = 


Qn 


; ^=[1 ij; x = lQn 0] 


(33) 


IV. Case Study 


This section validates the proposed generalized ESS model 
through time domain simulations. Three energy storage de¬ 
vices are considered, namely, SMES, CAES and BES. With 
this aim, the WSCC 9-bus test system (see Eig. Ell is used 
for all simulations. This benchmark network consists of three 
synchronous machines, three transformers, six transmission 
lines and three loads. Primary frequency and voltage regula¬ 
tors (AVRs) are also included. All dynamic data of the WSCC 
9-bus system as well as a detailed discussion of its transient 
behavior were originally provided in [30], and are publicly on 
several websites, e.g., [31]. 



Fig. 11: WSCC 9-bus test system with an ESS device connected to 
bus 8. 


An ESS is connected to bus 8. The capacities of the ESSs 
used in this case study have been designed according to the 
power installed in the system. This leads to assume large 
ESSs. 

Different scenarios have been performed in this paper: 
Subsection IIV- Al shows the response of a SMES connected to 
the WSCC system following a three-phase fault (IIV-A.II) . and 
stochastic variations of the loads (IIV-A.21 I. A similar analysis 
is carried out in Subsection IIV-BI for a CAES device. In this 
case, the contingency is a loss of load (IIV-B.II) . and stochastic 
processes are applied to all load power consumptions (IIV-B.2I) . 
Einally, a BES and a loss of load are considered in Subsec- 
tion llV-Cl Note that results shown in this section are not hand¬ 
picked but, rather, randomly selected among several hundreds 
of simulations that have been carried out to check the validity 
and accuracy of the proposed generalized ESS model. 

All simulations and plots have been obtained using DOME 
[32]. Dome has been compiled based on Python 2.7.5, 
CVXOPT 1.1.5, SuiteSparse 4.2.1, and Matplotlib 1.3.0; and 
has been executed on a 64-bit Linux Ubuntu 12.04 distribution 
running on 8 core 3.60 GHz Intel Xeon with 12 GB of RAM. 
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Time [s] 



Fig. 12: Response of the WSCC system with a SMES following 
a three-phase fault at bus 7. (a) Frequency of the COI. (b) Active 
power of the SMES. (c) Change in the Stored Energy of the SMES. 

A. SMES 

In this subsection, the storage device connected to bus 8 
is a 15 MW, 60 MJ SMES, whose model is described in 
Subsection IIII-BI To the best of our knowledge, the largest 
installed SMES (Center of Advanced Power Systems, Elorida 
State University) has a capacity of 100 MJ, and can provide 
100 MW peak and ±50 MW oscillatory power [10]. Two 
scenarios have been considered: Subsection IIV-A. l1 compares 
the dynamic response of the models of the SMES when the 



Time [s] 



Time [s] 



Time [s] 


Fig. 13: Response of the WSCC system following a three-phase fault 
at bus 7 when the SMES reaches its maximum storable energy, (a) 
Frequency of the COI. (b) Active power of the SMES. (c) Change 
in the Stored Energy of the SMES. 

system faces a three-phase fault, whereas Subsection IIV-A. 21 
considers stochastic variations of the loads and different initial 
states of charge of the SMES. 

1) Contingency: A three-phase fault occurs at bus 7 at t = 
1 s and is cleared after 70 ms through the disconnection of the 
line which connects buses 7 and 5. In this case, the frequency 
of the center of inertia (COI) of the system is regulated and 
its trajectory is shown in Eig. |12(a)| using the detailed and the 




























































































Fig. 14: Response of the WSCC system with a SMES considering stochastic variations of the loads, (a) The initial state of charge of the 
SMES is 20%. (b) The initial state of charge of the SMES is 50%. (c) The initial state of charge of the SMES is 80%. 


proposed models of the SMES. Figure [r2(a)| also shows the 
frequency of the COI when the SMES is modeled using the 
simplified model described in Subsection III-BI 

It can be observed that without ESS, the frequency variation 
after the fault is around 1%, and the steady-state is reached 
after about 40 s. The inclusion of the SMES in the system 
allows reducing the overshoot by 60%. Moreover, the settling 
time is about 15 s. 

Figure [T2(b)| compares the active power output of the SMES 
when the detailed and the proposed models are used. The 


SMES uses the load notation, therefore positive values of 
the power indicate that the SMES is storing energy, and vice 
versa. 

Finally, Fig. |12(c)| depicts the variation of the energy stored 
in the SMES. The base of the energy is 100 MJ, therefore the 
SMES increases its stored energy a maximum of about 57 MJ 
during and after the fault. The steady state value of the energy 
after the occurrence of the disturbance is about 20 MJ over 
the initial conditions. As it can be observed from Fig. [121 the 
performance of the system is basically the same when using 
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the detailed and the proposed models of the SMES. 

Figure [13] shows the performance of the system when 
the SMES reaches its maximum storable energy. In this 
simulation, it has been considered a maximum variation in the 
energy of 40 MJ over the initial condition. It can be observed 
that the effect of this sort of non-linearity is precisely captured 
by the proposed model, and the differences in the performance 
between the detailed and the proposed models are very small. 

Figures [T2| and [T3] also show the transient response of 
the SMES when the commonly-used simplified model of 
ESSs described in Subsection III-BI is applied. Gains and time 
constants of the different controllers defined in Sectionllllhave 
the same value for all ESS models. The trajectories obtained 
by imposing the simplified model consistently deviate from 
the behavior of the detailed one, particularly after the transient 
and after reaching the maximum variation of stored energy 
(note that, for a fair comparison, energy limits are also 
included in the simplified model in Fig. [T3T l. The proposed 
generalized model is able to accurately track the behavior 
of the detailed one in both, fast (transient) and slow (post¬ 
contingency) dynamics. On the other hand, the simplified 
model is accurate only in the first few seconds after the 
disturbance. The poor performance of the simplified model is 
an expected consequence of the reduced order of its dynamic 
equations. 

2) Stochastic Load Variations: In this scenario, all loads 
are modeled considering stochastic variations of the power 
consumptions. The Ornstein-Uhlenbecks process, also known 
as mean-reverting process, is used to model such varia¬ 
tions [33], [34]. The step size of the Wiener’s process is 
h — 0.01s, and the time step of the time integration At is 
set to 0.1s. Initial values of the load and generation are set 
using a uniform distribution with 5% of variation with respect 
to the base-case. 

Figure [T4| shows the frequency of the COI and the state 
of charge of the SMES for three different load profiles and 
for an initial state of charge of 20%, 50%, and 80% in 
Figs. |14(a)| |14(b)[ and |14(c)l respectively. These percentages 
are expressed in terms of the allowed energy variability of 
the storage device. For each simulation, SMES models are 
simulated considering same load variation profiles. 

The following remarks are relevant: 

i. The proposed generalized ESS model tracks the beha¬ 
vior of the detailed one better than the commonly-used 
simplified model of the SMES, for any initial state of 
charge. 

ii. The proposed model is more accurate the closer is 
the state of charge to the initial condition during the 
simulations. This is an expected property of linearized 
models. 

iii. Based on the results of hundreds of simulations with 
different load profiles and initial conditions, it has been 
observed that the best average accuracy of the proposed 
generalized model is obtained by computing the matrices 
of ([T2l) - ([T3l) for an initial state of charge of 50%. 

B. CAES 

In this case study, we consider a 15 MW, 30 bar CAES, 
whose model is described in Subsection IIII-CI The CAES 
regulates the active power flowing through the transformer 
connecting buses 2 and 7. As in previous Subsection, 



Time fsl 



Time [s] 

Fig. 15: Response of the WSCC system with a CAES following a 
loss of load, (a) Power flowing through the transformer connecting 
buses 2 and 7. (b) Active power of the CAES. 

the CAES uses the load notation, and two scenarios are 
performed: Subsection IIV-B.II considers a loss of load as 
contingency, while Subsection IIV-B.2I includes stochastic 
variations in all loads. Note that the maximum size of 
currently installed above-ground CAESs is up to 5 MW [1]. 
In this example, we assume that the CAES has three times 
this capacity. 

1) Contingency: The contingency is a loss of load of 15 
MW occurring at bus 5 at f = 10 s. Finally, the load is 
reconnected at f = 80 s. Figures |15(a)| and |15(b)| illustrate 
the active power flowing through the transformer connecting 
buses 2 and 7 and the active power of the CAES, respectively. 
It can be seen that the response of the CAES simulated using 
the proposed model is accurate despite the complexity of the 
detailed CAES model shown in Subsection IIII-CI 

2) Stochastic Load Variations: In this scenario, all loads 
are modeled using the mean-reverting process, (see Subsec¬ 
tion IIV-A. 2 I 1 . Figure |16(a)| depicts the active power of the 
CAES for each model of this device, for an initial state of 
charge of about 50%. All control parameters are the same for 
all models. The time constants of the simplified model of the 
CAES have been tuned after several trial and error attempts, in 
order to obtain the behavior as close as possible to the detailed 
model. Figure [T6(a)| shows that the response of the commonly- 
used simplified model differs considerably from the one of the 
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Fig. 16: Response of the WSCC system with a CAES considering 
stochastic perturbations in the loads, (a) Same control parameters for 
all models, (b) Control parameters of the simplified CAES model are 
three times smaller than those of the detailed and the generalized 
ones. 

detailed model. On the other hand, the proposed generalized 
model appears to be very accurate. In order to obtain a more 
accurate response of the storage device using the simplified 
model, it is required to properly tune the control gains Kpp 
and Kip of the controller depicted in Fig. |7] Figure |16(b)| 
depicts the response of the CAES considering the same load 
profiles as in Fig. |16(a)| and when the gains ATpp and Kip 
are three times smaller than the gains Kp^ and Ki^ of the 
storage controller of Fig. |6] respectively. Some remarks can 
be deduced: 

i. The accuracy of the proposed generalized model does not 
appear to be affected by the complexity and the size of 
order reduction of the original detailed model. Note that 
the order reduction achieved in the case of the CAES is 
from 29 variables in the detailed model to only 5 in the 
proposed one. 

ii. The accuracy of the commonly-used simplified model 
cannot be guaranteed even if a lengthy and careful tuning 
of its control parameters is carried out. 

iii. Because of the required tuning, the design of control 
strategies cannot be based on the simplified model. On 
the other hand, exactly same control parameters can be 
used for both the generalized and detailed ESS models. 



Time [si 



Time [s] 


Fig. 17: Response of the WSCC system with a BES following a loss 
of load, (a) Erequency of the COI. (b) Active power of the BES. 


C. BES 

In this example, a 40 MW BES is installed at bus 8 [13]. 
The data of the BES is taken from [35], that describes 
a 55 MW, 76.7 GJ BES used for frequency control. The 
contingency is a loss of a 40 MW load at bus 5. The load is 
disconnected at f = 10 s, and is reconnected at f = 100 s. In 
this case study, the BES regulates the frequency of the COI. 

Eigures |17(a)| and |17(b)| show the frequency of the COI 
and the power output of the BES, respectively. The initial 
SOC of the battery is set to 85% to force operating close to 
the nonlinear voltage-current characteristic of the battery [12], 
[13]. As stated in Subsection IIII-DI two sets of matrices for 
the general model of the BES have to be considered because 
of the discontinuity in the polarization voltage, Up, depending 
on whether the battery is charging (time from f = 10 s to 
t = 100 s) or discharging (time from t = 100 s). It can be 
seen from Eig. [17] that the generalized model is able to track, 
with a good level of accuracy, the behavior of the detailed 
model of the BES for both operating conditions, despite the 
nonlinearities of the device at such a high initial SOC, and 
the switching between the two operating modes. 

D. Concluding Remarks 

The following remarks on the proposed generalized model 
of ESSs are relevant. 
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1) The proposed model provides a good compromise bet¬ 
ween simplicity and accuracy. In fact, the proposed 
model has a reduced and constant dynamic order but, 
nevertheless it can reproduce faithfully the behavior 
of ESSs whose detailed transient stability models are 
considerably more complex (e.g., CAES). On the con¬ 
trary, the case study shows that simpler models of ESSs 
might not be precise, particularly if windup limiters are 
binding. 

2) Linearization of a subset of ESS equations does not 
affect accuracy. The rationale behind this observation is 
that most ESS variables involved in nonlinear equations 
are bounded and, hence, the ESS operating point does 
not vary consistently even after a large perturbation. 

3) Detailed ESS models are needed to define the para¬ 
meters of the proposed model. These data have to be 
provided by ESS manufactures, in the same way makers 
provide d- and q-axis reactances and time constants of 
the synchronous machine Park model. Hence, we con¬ 
sider that the proposed ESS model can be an opportunity 
to define a ESS standard for transient stability analysis. 

4) Control strategies designed for the proposed generalized 
model can be directly applied to the detailed ESS 
models. The linear structure of the equations of the 
proposed model simplifies the design of more advanced 
and robust control strategies for ESSs that could be 
applied subsequently to the detailed models of these 
devices. 

V. Conclusions 

This paper proposes a generalized model of energy storage 
devices for voltage and angle stability analysis. Such a model 
allows simulating different technologies using a fixed set 
of DAEs and parameters. Simulation results show that the 
proposed model is able to accurately reproduce the dynamic 
behavior of detailed transient stability models for large dis¬ 
turbances, namely faults and loss of loads, and across their 
whole operating cycle. The non-linearity of ESS controllers, 
i.e., hard limits, are also properly taken into account by the 
proposed model. 

The proposed generalized ESS model appears to be par¬ 
ticularly useful to synthesize and compare different control 
strategies. In fact, since the proposed model is composed of 
a fixed set of equations, the same control scheme can be 
straightforwardly used for testing the dynamic response of 
different technologies. This will be addressed in future work. 
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